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Abstract 



We introduce a new class of non-isothermal models describing the evolution 
of nematic liquid crystals and prove their consistency with the fundamental 
laws of classical Thermodynamics. The resulting system of equations captures 
^ ■ all essential features of physically relevant models, in particular, the effect of 

stretching of the director field is taken into account. In addition, the associated 
I initial-boundary value problem admits global-in-time weak solutions without 

any essential restrictions on the size of the initial data. 

o 

1 Introduction 

The celebrated Leslie- Ericksen model of liquid crystals, introduced by Ericksen [6] and 
^ ' Leslie [H], is a system of partial differential equations coupling the Navier-Stokes equa- 

■ tions governing the time evolution of the fluid velocity u = u(t, x) with a Ginzburg- 

Landau type equation describing the motion of the director field d = d{t, x) , repre- 
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senting preferred orientation of molecules in a neighborhood of any point of a reference 
domain. 

A considerably simplified version of the Leslie- Ericksen model was proposed by 
Lin and Liu [15], [16], and subsequently analyzed by many authors, see [20], [21], [25] 
among others. The simplified model ignores completely the stretching and rotation 
effects of the director field induced by the straining of the fluid, which can be viewed 
as a serious violation of the underlying physical principles. 

Such a stretching term was subsequently treated by Coutand and ShkoUer [1], 
who proved a local well-posedness result for the corresponding model without thermal 
effects. The main peculiarity of this model is that the presence of the stretching term 
causes the loss of the the total energy balance, which, indeed, ceases to hold. In order 
to prevent this failure. Sun and Liu [23] introduced a variant of the model proposed by 
Lin and Liu, where the stretching term is included in the system and a new component 
added to the stress tensor in order to save the total energy balance. A more general 
class of models based on the so-called Q-tensor formulation was recently introduced 
in [21 [18] in the isothermal case. 

Motivated by these considerations, in the the present contribution, we propose a 
new approach to the modeling of non-isothermal liquid crystals, based on the principles 
of classical Thermodynamics and accounting for stretching and rotation effects of the 
director field. To this end, we incorporate the dependence on temperature into the 
model, obtaining a complete energetically closed system, where the total energy is 
conserved, while the entropy is being produced as the system evolves in time. We 
apply here the mechanical methodology of [TU], which basically consists in deriving 
the equations of the model by means of a generalized variational principle. This 
states that the free energy \1/ of the system, depending on the proper state variables, 
tends to decrease in a way that is prescribed by the expression of a second functional, 
called pseudopotential of dissipation, that depends (in a convex way) on a set of 
dissipative variables. In this approach, the stress tensor a , the density of energy vector 
B and the energy flux tensor H are decoupled into their non-dissipative and dissipative 
components, whose precise form is prescribed by proper constitutive equations (see 
below for details). It is interesting to note that the form of the extra stress in the 
Navier-Stokes system obtained by this method coincides with the formula derived 
from different principles by Sun and Liu in [23]. 

The system of partial differential equations resulting from this approach couples 
the incompressible Navier-Stokes system for the velocity u, with a Ginzburg-Landau 
type equation for the director field d and a total energy balance together with an en- 
tropy inequality, governing the dynamics of the absolute temperature 6 of the system. 

Leaving to the next Section [2] the complete derivation of the model, let us just 
briefly introduce here the PDE system we deal with. The Navier-Stokes system couples 
the incompressiblity condition 

divu = (1.1) 

with the conservation of momentum 

Ui u ■ V^n + VxP = divS + divcr"'^ + g, (1.2) 
where p is the pressure, and the stress is decomposed in a dissipative and non dissi- 
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pative part, respectively given by 



S = ^ (V,u + V*.u) , a"'^ = -AV,d V,cZ + A(/(d) - Ad) ® d, (1.3) 

where we have set Vxd Vxd := J2k didkdjdk ■ 
The director field equation has the form 

dt + uVxd-d- V,u = 7(Ad - /(d)), (1.4) 

where /(d) = ddF{d) and F penalizes the deviation of the length \d\ from the value 
1. It is a quite general function of d that can be written as a sum of a convex (possibly 
non smooth) part, and a smooth, but possibly non-convex one. A typical example is 
F(d) = (|dp-1)2. 

Finally, the total energy balance 



Q|up + + u ■ Va; + + div (^pu + g - §u - a"'^uj (1.5) 



= ^ ■ u + A7div {Wxd ■ {Ad - f{d)) ) , 
with the internal energy and the flux 

e = ^|V^d|^ + AF(d) + ^, q = q'^-AV^d-V^u-d, q'^ = -k{e)\/J-h{e){d-\/J)d, 

is coupled with the entropy inequality 

H{9)t + u ■ VxH{9) + div {H'{9)q'^) (1.6) 

>H'{e) {S:VxU + X-f\Ad-f{d)\^)+H"{e)q''-VJ, 

holding true for any smooth non-decreasing concave function H . The derivation of 
the above system will be detailed in the next Section [21 while the remainder of the 
paper will be devoted to the proof of a global existence result for the corresponding 
initial-boundary value problem in the framework of weak solutions in f2 x (0, T) , being 
Q a bounded and sufficiently regular subset of M"^ and T a given final time. 

Let us note that the model obtained here looks quite different from the one 
obtained in [5]. This is mainly due to the presence in the internal energy e of the 
quadratic term |Va;dp (that is related to the expression (12.11) of the free energy func- 
tional) and to the stretching term d ■ V^u in (11.41) which produces, in order that the 
principles of Thermodynamics are respected, two new non dissipative contributions in 
the stress tensor § in (11.31) and in the flux q . Actually, the latter is given here by the 
sum of a standard heat flux and of an elastic part given by the term — AV^d ■ V^u ■ d 
(cf. the next Section [2] for further details on this point). 

Indeed, in contrast with [9] , the presence of the stretching term d ■ V^^u in the 
director fleld equation prevent us from applying any form of the maximum principle to 
(11.41) . Hence, we cannot recover an L°° -bound on d (which we obtained, instead, in 
[H]). However, we can still get here the global existence of weak solutions to the initial 
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boundary value problem coming from the PDE system fll.lH1.6p without imposing any 
restriction on the space dimension, on the size of the initial data or on the viscosity 
coefficient /i (such a restriction was taken in the paper [23], devoted to an isothermal 
model closely related to ours). In this sense, our results can be seen as a generalization 
of those obtained in [23]. 

The compatibility of the model with First and Second laws of thermodynamics 
turns out to be the main source of a priori bounds that can be used, in combination 
with compactness arguments, to ensure stability of the family of approximate solutions. 
The key point of this approach is replacing the heat equation, commonly used in models 
of heat conducting fluids, by the total energy balance (11.51) . Accordingly, the resulting 
system of equations is free of dissipative terms that are difficult to handle, due to the 
low regularity of the weak solutions. In contrast with the standard theory of Navier- 
Stokes equations, however, we have to control the pressure appearing explicitly in the 
total energy flux and, in order to do that we will need to assume the complete slip 
boundary conditions on the velocity flled u (cf. (13. ip ). Note that a similar method 
applied to different models has been recently used in [1], [7], and [9]. 

Finally, let us notice that the non-isothermal liquid crystal model accounting for 
the stretching contribution has also recently been analyzed in [3] (in case of Dirichlet 
boundary conditions for u and Neumann or non homogeneous Dirichlet boundary 
conditions for d) and in [19], where the long time behaviour of solutions is investigated 
in two cases: in the 3D case without any condition on the size of the viscosity coefficient 
H and in case of a non analytic nonlinearity / . Both these results generalize the ones 
obtained in [25] . 

The paper is organized as follows. In Section [21 we give a detailed derivation 
of the model and discuss its compatibility with the basic laws of Thermodynamics. 
In Section [31 we introduce some technical hypotheses and formulate the main result 
concerning existence of global-in-time weak solutions to the resulting PDE system. 
Finally, the last two Sections [4] and [5] are devoted to the proof of the existence result 
via approximation, a-priori estimates and passage to the limit techniques based on 
lower semicontinuity and convexity arguments. As already pointed out, the energy 
balance is written in the form of a conservation law for the total energy rather than 
for the temperature, where the highly non-linear terms dissipative terms are absent. 
The price to pay is the explicit appearance of the pressure in the global energy balance 
determined implicitely by the Navier-Stokes system. 



2 Mathematical model 

We suppose that the fluid occupies a bounded spatial domain 17 C , with a suffi- 
ciently regular boundary, and denote by u = u(t, x) the associated velocity field in the 
Eulerian reference system. Moreover, we introduce the absolute temperature 6(t,x) 
and the director field d{t, x) , representing the preferred orientation of molecules in a 
neighborhood of any point of the reference domain. Furthermore, we denote 

dw 

— = W = Wt + U- \/^w, 
dt 



the material derivative of a generic function w, while Wt (or also dfw) denotes the 
partial derivative with respect to t. 
Finally, the quantity 

^ = dt + \l- Vxd - d ■ Va:U 

characterizes the total transport of the orientation vector d. Note that the last term 
accounts for stretching of the director field induced by the straining of the fluid. 



2.1 Free-energy and pseudopotential of dissipation 

Following the general approach proposed in the monograph [TU] , we start by specifying, 
in agreement with the principles of classical Thermodynamics, the free-energy and the 
pseudopotential of dissipation. The interested reader may consult [TOl Chapters 2,3] 
for details. 

We begin by introducing the set of state variables, describing the actual config- 
uration of the material, specifically, 

E = {d,v.,d,e). 

Next, the set of the dissipative variables describing the evolution of the system, and, 
in particular, the way it dissipates energy, is given by 

where 

(V.U + Vi.u) 



2 

denotes the symmetric gradient of u. 

Motivated by the original (isothermal) theory proposed by Ericksen [5] and 
Leslie [13], we choose the free energy functional in the form 

^{E) = ^|v.d|2 + XF{d) -e\oge, (2.1) 

where A is a positive constant. The function F in (12. ip penalizes the deviation 
of the length \d\ from its natural value 1; generally, F is assumed to be a sum 
of a dominating convex (and possibly non smooth) part and a smooth non-convex 
perturbation of controlled growth. A typical example is F{d) = (jcZp — 1)^. For the 
sake of simplicity, we have assumed that the thermal and elastic effects are uncoupled 
in 

The evolution of the system is characterized by a second functional $, called 
pseudopotential of dissipation, assumed to be nonnegative and convex with respect to 
the dissipative variables. Specifically, we consider $ in the form 



HSE, E) = ^|.(u)r + /o(divu) + ^|V.^r + I 



Dd 



Dt 



20 \ ^1 
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where fi = fi{9) > is the viscosity coefficient, 77 > is a constant, and k, h 
represent the heat conductivity coefficients - positive functions of the temperature. 
The incompressibility of the fluid is formally enforced by /q - the indicator function of 
{0} (given by Jq = if divu = and +00 otherwise). 

2.2 Constitutive relations 

We start by introducing the stress tensor a , the density of energy vector B , and the 
energy flux tensor H; all assumed to be the sum of their non-dissipative and dissipative 
components, namely, cr = cr"'^ + a"^ , B = B"*^ + B*^ , H = M^'^ + , where 

. d-^ dF 

B'=^ = (2^3) 

' Dt' ^ ' 

M"'^ = — — = XV,d. (2.4) 

Moreover, we set H*^ = 0. 

The heat and entropy fluxes (denoted respectively by q"' and Q) are 

q'^ = = -9^^ = -mV^e - h{e){d ■ VJ)d. (2.5) 

The stress tensor a consists of two parts: the dissipative one 

a''=^^ = ^i{d)e{u)-pl=:^-pl, (2.6) 

-p G 9/o(divu), § = ^{9)e{\i), 

and the non dissipative part a"'^ to be determined below (cf. f l2.10p and f l2.12p ). 
The entropy of the system is given by 

s = -^ = l + \oge (2.7) 

and, finally, the internal energy e reads 

e=^ + es = e + XF{d) + ^\S/^d\'^. (2.8) 

2.3 Field equations 

In accordance with Newton's second law, the balance of momentum reads 

9tu + div(u(g)u) = diver + sf, (2.9) 
where gr is a given external force. 
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The entropy balance can be written in the form 

st + n- V,.s + div Q = i (^a'' : e{u) + B''-^-Q- VJ^ (2.10) 
or equivalent ly 

+ divq'^ = : e{u) + B'^ ■ (2.11) 

In agreement with Second law of Thermodynamics, the right hand side of (12.101) is 
non-negative. 

The balance of internal energy reads 

^ n , ^ ^ Dd „ Dd , , 

+ u ■ V^e + divq = a : e(u) + B ■ — + H : V^— , (2.12) 

with the internal energy flux q = q'^ + q^'^ , where the dissipative part q'^ is given by 
(12. Sp . while the non-dissipative component will be determined below. 

Finally, the equation which rules the evolution of the orientation vector d is 
derived from the principle of virtual powers (cf. [101 Chap. 2]) and it takes the form 

divH-B = 0, (2.13) 

specifically, 

+ u ■ V.d - d ■ V.u = 7(Ad - /(d)), 7 = X/'q. (2.14) 

The non-dissipative component of the stress cr"^ and of the flux q^'^ are deter- 
mined by means of (12.101) . (I2.12p . and the constitutive relations derived above. Indeed, 
computing ^ by means of the standard Helmholtz relations, we get 



whereas 



de ^ds dO ^ dd ^ diVxd) ^ds . 



Thus, rewriting fl2.15p with help of fl2.16p . and expressing 6^ by means of (12. lip , we 
get, thanks also to ([22H2S]), 

divg'^'^ - a"" : V.u = -A(/(d) ® d) : V.u - \Y,d.Md%,.,u,)d,. (2.17) 
Therefore, 

a'^'^ = -AV^d V^.d + A(/(d) - Ad) ® d, q'"^ = -AV^d ■ V^u • d. (2.18) 

Summing up the previous discussion, we arrive at the following system of equa- 
tions: 

INCOMPRESSIBILITY: 

divu = 0; (2.19) 
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CONSERVATION OF MOMENTUM: 
uj + u ■ + VxP = div S + div cr"'^ + (2.20) 
where p is the pressure, and 

S = ^ (V,u + V^u) , a"'^ = -AV,cZ © + A(/(d) - Ad) ® d; (2.21) 

DIRECTOR FIELD EQUATION: 

dt + u- V^d - d ■ V^.u = 7(Ad - /(d)), (2.22) 
TOTAL ENERCY BALANCE: 



dt Q|u|2 + e^ +u- V,; Q|up + e^ + div (^pu + q - Su - a""'u) 
= ^ ■ u + A7div (V,d ■ (Ad - /(d)) 



(2.23) 



with the internal energy 

e = ^|V^d|2 + AF(d) + e 

and the flux 

q = q-^ + q"'^ = -k{9)VJ - h{9){d ■ VJ)d - AV^.d ■ V^^u ■ d, 
together with 

ENTROPY INEQUALITY: 

H{e)t + u ■ V,H{9) + div {H'{d)q'') (2.24) 

> H'{9) (S : V.u + A7I Ad - /(d) p) + H"{9)q'' ■ VJ, 

holding for any smooth, non-decreasing and concave function H . 

Actually, the total energy balance 02.231) follows easily from (12. 9p . f l2.12p . com- 
bined with (12.131) . It is remarkable that equations ([2IISH222D in the isothermal case 
reduce to the model derived by Sun and Liu in [23] by means of a different method. 

3 Main results 

3.1 Initial and boundary conditions 

In view of a rigorous mathematical study, system (l2.19f[T^2]) must be supplemented 
by suitable boundary conditions. Actually, to avoid the effect of boundary layer on 
the motion, we assume complete slip boundary conditions for the velocity: 

u • n\9n = 0, [(S + a"^)n] x n|an = 0. (3.1) 
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Moreover, we consider no-flux boundary condition for the temperature 

■ nlan = 0, (3.2) 

and Neumann boundary condition for the director field 

Vxdi ■ nlan = for i = 1,2,3. (3.3) 

The last relation accounts for the fact that there is no contribution to the surface 
force from the director d. Note that the above conditions are also suitable for the 
implementation of a numerical scheme (see [17] for further comments on this point). 
Of course, we also need to assume the initial conditions 

u{0, ■) = Uo, d{0, ■) = do, 9{0, ■) = ^0, (3.4) 

In the remaining part of the paper, our aim will be that of showing existence of global- 
in-time solutions to system fl2.19f[T^^ . coupled with the above initial and boundary 
conditions and without assuming any essential restriction on the data. 



3.2 Weak formulation 

In the weak formulation, the momentum equation (12.201) . together with the incom- 
pressibility constraint fl2.19p . and the boundary conditions (13.11) . are replaced by a 
family of integral identities 

/ u{t, ■) ■ V^^^ = for a.a. t G (0, T) (3.5) 
Jn 

for any test function (f G C°°{Q), and 

/ / ( ■ + uC^ u : + p div tp 
Jo Jn ^ 



nd\ 



(3.6) 



5 + a"''):V,^- / g-^ - / tio-<^(0,-), 
Jo. Jn Jn 

for any Lp E C^{[0, T) xTt;^?) , if- n\Qn = 0. Note that (^M) includes also the initial 
condition u(0, ■) = uq. 

Equation fl2.22p describing the evolution of the director field d will be satisfied 
in the strong sense, more specifically, 

dtd + u-V^d- d-VxU = -f(^Ad- f{d)^ a.e. in (0, T) x fi, (3.7) 

together with 

V xdi ■ nidn = 0, i = 1,2, 3, d(0, ■) = do- 
Similarly, the weak formulation of the total energy balance f l2.23|) reads 



+ [ [ {pu + q-Su- (t""'u) ■ W^ip 
Jo Jn 

= A7^ ^(v,.d-(Ad-/(d))).V,.^-^ ^^■uv9-^Q|nop + eo)<^(0 

for any (f G C^{[0, T) xQ), where eo = ||V^do|^ + AF((io) + ^o- 
Finally, the entropy inequality f l2.24p is replaced by 



r [ H{9)dt^+ f [ {H{d)u + H\d)q'')-V^^ (3.9) 
JO Jn Jo Jn 

f [ {H\e) (§: V.u + A7|Ad-/(d)|2)+ff"(^)g^.V.^)(^- [ HiOoMO,-) 
Jo Jn Jn 

for any if G C^([0,T) x Q) , ip>0, and for any smooth, non-decreasing and concave 
function H. 

A weak solution is a triple (u, d, 9) satisfying f l3.5H3.9p . 
3.3 Main existence theorem 

Before formulating the main result of this paper, we list the hypotheses imposed on 
the constitutive functions. Specifically, we assume that 

F e C^(M^), F > 0, F convex for all \d\ > Dq, lim F{d) = oo, (3.10) 

for a certain Dq > 0. 

The transport coefficients fi, k, and h are continuously differentiable functions 
of the absolute temperature satisfying 

< /i < M < /J, < ^ < k{9), h{9) < k for all ^ > (3.11) 

for suitable constants k, k, /i, JI. 

Our main result reads as follows. 

Theorem 3.1. Let C M'^ be a bounded domain of class C'^^'^ for some u > 0, 
g G L^((0,T) X n;R^). Assume that hypotheses (E2S|), (EUP are satisfied. Finally, 
let the initial data be such that 

uo G L2(fi;M3)^ div^^Q = 0, do G W^'\n;R^), F{do) G L\n), 

(3.12) 

9o G L^i^l), essinfQ6'o > 0. 

Then, problem (\3.5[]3.9\] possesses a weak solution (u, d, 9) in (0, T) x VL be- 
longing to the class 

u G L~(0, T- L\n; R^)) n ^^(0, T; W^'^Q; R^)), (3.13) 

d G L°°(0, T; W^'\n; R^)) n ^^(0, T; W^'^Q; R^)), (3.14) 



< 
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F{d) G L°°(0, T; L\n)) f] 1^/^(0, T) x fi), (3.15) 

e G L°^{0,T;L\n))nLP{0,T;W^'P{n)), l<p< 5/4, ^ > a.e. in (0,T)xfi, (3.16) 
with the pressure p , 

pG L^/^((0,T) X fi). (3.17) 
The rest of the paper is devoted to the proof of Theorem 13.11 

4 A priori bounds 

In this section, we collect the available a priori estimates. These will assume a rigorous 
character in the framework of the approximation scheme presented in Section [S] below. 

Integrating (12.231) over Q and using Gronwall's lemma, we immediately obtain 
the following bounds: 

u e L'^{0,T;L\n;R^)), (4.1) 

9 e L°^{0,T;L\n)), (4.2) 

d G L°°(0, T; W^'\Q; R^)), F{d) G L°°(0, T; L\Q)), (4.3) 

where we have used hypotheses fl3.10p . (13.111) . 

Similarly, integrating (I2.24p with H{6) = 6, and using (14. 2p . we obtain 

e{u) G L2((0, T) X Q, M3x3)^ Ad - /(d) G L^{{0, T) x VL; M^). (4.4) 

yielding, by virtue of (14. ip and Korn's inequality, 

u G L2(0, T; W^^'^iVt- R?)) n L^°/^((0, T) x fi; M^). (4.5) 

Moreover, it follows from (14.41) and convexity of F (cf. hypothesis (13.101) ) that 

/(d)GL2((0,T)xf];R3); (4.6) 

therefore, using (14. 4p again we infer that 

d G L2(0,T;iy^'2(n;M3))_ (4 7) 

Interpolating (14. 3 p and (14. 7p we get 

d G L^°((0, T) X fi; M=^), V^d G L^°/^((0, T) x fi; M^xS)^ 
whence (cf. ([221])) 

By the same token, by means of convexity of F (cf. (13.101) ). we have 

|F(d)|<c(l + |/(d)||d|), 

yielding 

F(d) G L5/3((0,T) X (]). (4.9) 
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As the velocity satisfies the shp boundary conditions fl3.ip . the pressure p can 
be "computed" directly from (12.201) as the unique solution of the elliptic problem 

Ap = div div (s + a'"^ - u (g) u j + div g, 

supplemented with the boundary condition 

VxP ■ n = (div [s + a^"^ - u^u) + g) ■ n on dQ . 

To be more precise, the last two relations have to be interpreted in a "very weak" 
sense. Namely, the pressure p is determined through a family of integral identities: 

/ pA<^ = / (§ + d"'^ - -u ® : Vl^ - g- V^y?, (4.10) 
Jn Jn ^ ' Jn 

for any test function ip G C°°{Q), V^v^ ■ nlgQ = 0. Consequently, the bounds estab- 
lished in (14. 5 p and (14. 8 p may be used, together with the standard elliptic regularity 
results, to conclude that 

p G L5/^((0,T) X fi). (4.11) 

Finally, the choice H{9) = (1 + 9^, rj e (0,1), in ([2Jl]), together with the 
uniform bounds obtained in (l4.1H4.5p . yields 

V^(l + ey e L\{0, T) X Q; M^) for any < < ^ . (4.12) 

Now, we apply an interpolation argument already exploited in [T]. Using (14. 2 p and 
fl4:T2|) and interpolating between 6 G L'^iO^T; L\n)) and 6" e L^(0, T; L3(fi)) , for 
I' G (0, 1) , we immediately get 

9 G L'?((0, T) X n) for any 1 < g < 5/3 . (4.13) 

Furthermore, seeing that 



J{o,T)xn \J{o,T)xn J \J{o,T)xn 

for all p G [1, 5/4) and z^ > 0, we conclude from (I4.12p and (14.130 that 

G LP((0,T) X n;R^) for any 1 < p < 5/4. (4.14) 
Finally, the same argument and H{9) = log 6* in (I2.24p give rise to 

log^ G L\{0,T)-W^'\n)) n L'=^{0,T- L\n)), (4.15) 

where we have used (14. 2p . 

The a priori estimates derived in this section comply with the regularity class 
(l3.13f[XT7|) . Moreover, it can be shown that the solution set of (l3.5H3.9p is weakly sta- 
ble (compact) with respect to these bounds, namely, any sequence of (weak) solutions 
that satisfies the uniform bounds established above has a subsequence that converges 
to some limit that still solves the system. Leaving the proof of weak sequential stabil- 
ity to the interested reader, we pass directly to the proof of Theorem 13.11 constructing 
a suitable family of approximate problems. 
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5 Approximations 

For the sake of simplicity, we restrict ourself to the case g = and A = 7 = 1. 
Solutions to the Navier-Stokes system (13. 5p . fl3.6p will be constructed by means of 
the nowadays standard Faedo-Galerkin approximation scheme, see Temam [23]. Let 
W^'^(f2; M^) be the Sobolev space of solenoidal functions satisfying the impermeability 
boundary condition, specifically, 

W^ll = {v E W^'^in; R^) I divi; = a.e. in fi, v ■ n\an = 0} . 

Since dQ is of class C'^'^'^ , there exists an orthonormal basis {f^j^^i of the Hilbert 
space W^'l such that G C^^" , see Theorem 10.13]. We take M < and 
denote = spa.n{vn}n=i, and [v]m - the orthogonal projection onto the space 

The approximate velocity fields wat^a/ G C^([0, T]; Xat) solve the Faedo-Galer- 
kin system 



dt 



/ Un,m-v= [ [un^m]m<^Un,m -"^xV--^ [ \V^un,m\'' '^^xUnm - ^xV (5.1) 
Jo, Jn ^'^ Jn 

- J ^^^^'^^^ (^'^xUn,m + V^un^m^ : V^t; + j VxdNM © ^xdNM ■ ^xV 



- I {f{dN,M) - ^dN,M) ® dN,M ■ ^xV, 



n 

/ *t7V,M(0, ■)-V = / Uo-V, 

Jn Jn 

for any v G X^, where r G (3,10/3). The extra term j^lV xUn,m\^~'^'V xUn,m guar- 
antees sufficient regularity for the velocity field needed in the director equation. Our 
strategy is to pass to the limit first for N ^ 00 and then for M — >■ 00 . 

The functions dj^^M are determined in terms of Uat^m as the unique solution of 
the parabolic system 

dtdN,M + Un,M ■ '^xdN,M — dN,M ' '^xUn,M = '^d^^M — fid^^M), (5.2) 

supplemented with 

V^.(djv,M)i ■ n\dn = 0, z = 1, 2, 3, (5.3) 
d7V,Af(0, ■) = c«o,M, (5.4) 

where do^M is a suitable smooth approximation of do . 

Next, given un,m, d^^M , the temperature 9n,m is determined as the unique 
solution to the heat equation (cf. Ladyzhenskaya et al. [TTl Chapter V, Theorem 8.1]): 

dtON,M + div {9n,mUn,m) + div q^^^ (5.5) 
1 ^ 2 

Qif,M ■ ^\dn = 0, (5.6) 
13 



^^,Af(0,-) = (5-7) 



where E>n,m = ^i^^f^ ( V^uat^m + V^uat^m) , and 

Actually, relation (15.51) is just an explicit reformulation of (12.111) . 

Finally, the pressure pn,m is found as before as the (unique) solution to a system 
of integral identities: 

Pn,m^'^ (5-8) 



(§Af,M-Va;dAr,M©Va;dAr,M-(AdAr,M-/(t^Af>/))®t^Af>/)-[nAr,M]M®'"iV,M) : V^^J 
+ ^ y" \'^xUn,m\''~'^'^xUn,M ■ V^v^, 

satisfied for any test function ip G C°°{Q), Vxf ■ n\Qn = 0. 

Regularizing the convective terms in (15. ip is in the spirit of Leray's original 
approach [12] to the Navier-Stokes system. As a result, we recover the internal energy 
equality at the level of the limit — oo. This fact, in turn, enables us to replace 
the internal energy equation (15.51) by the total energy balance before performing the 
limit M — )■ oo . For fixed M, A^ , problem (I5.1H5.8P can be solved by means of a 
simple fixed point argument, exactly as in [H Chapter 3]. Note that all the a priori 
bounds derived formally in Section H] apply to our approximate problem. Thus, given 
u e C([0,T]; ATjv), we can find d = d[u] solving ([S2HS3]), and then 9 = e[u,d] 
and the pressure p satisfying (l5.5H5.8p . Plugging these functions d, ^ in (15. ip . the 
corresponding solution T[u] then defines a mapping u i— t- Tlu] . By the a priori 
bounds obtained in Section |U we can easily show that T possesses a fixed point by 
means of the classical Schauder's argument, at least on a possibly short time interval. 
However, using once more the a priori estimates we easily conclude that the approxi- 
mate solutions can be extended to any fixed time interval [0, T] , see [HI Chapter 6] for 
details. 



5.1 Passage to the limit as N ^ oo 

Having constructed the approximate solutions un,m , d^^M , 6n,m^ ^^"^ Pn,m, we let 
N ^ oo. To take the limit, we need to modify a bit the formal estimates obtained in 
Section m taking care of the regularizing terms added in (15. ip and (15. 5p . Indeed, from 
the energy estimate we now additionally obtain 

M ^||Va;'UAr,M||Lr((o,T)xn;R3x3) < C , (5.9) 

whence we infer that {V xUn^m\'^~'^'^ x'i^N,M is uniformly bounded in L~ {{0,T) x Q) 
for fixed M . Moreover, in place of (14.111) we deduce from (15.81) the estimate 

|b7V,M||L^/-i((0,T)xf7) < C{M) , (5.10) 
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where we observe that 

r /lO 3\ / 10\ , , 

— 5). sincere (3.-). (5.11) 

Note that, at least at the level of approximate solutions, relation f l2.24p holds 
true as an equality. Hence, taking H{9) = (1 + O)"^ , with r] G (0, 1) , in f l2.24p . we get 

ll^t^Af,Afll{CO([0,T];Wl.''(Q)))* < C\\dt6'^^]^j\\Li({0,T)xn) < C, 

where C is a positive constant independent of and M, with s G (3, +oo) , u G 
(0, 1/2) . This leads to the convergence relations: 

Un,m Um weakly-(*) in L°°(0, T; L'^{n; M^)) n ^^(0, T; W^'^Q; R^)) , (5.12) 
V^U7v,M ^ Vx.um weakly in L%0, T; U{VL; M^)) , (5.13) 

dtVLNM dtVLM weakly in L2(0,T; (iy^'2(fi; M^))*) + L^(0,T; PF-^''^/"-^(fi; M^)) , 

(5.14) 

Pn,m ^ Pm weakly in L^/^^i((0, T) x fi) , (5.15) 
^^,Af ^ weakly-(*) in ^^(0, T; W^^^iSl)) n L°°(0, T; Li/'^(fi)) , (5.16) 
dte^^M -> 5*^A/ weakly-(*) in (Co(0, T; lyi'^(fi)))* , (5.17) 
log ^iv,M ^ log^Af weakly in L2(0,T; ^^^'^(fi)) , (5.18) 
dN,M ^ dA/ weakly-(*) in L°°(0,T; iy^'2(fi;M3)) p L2(0,T; W'^''^{yt-W^)) , (5.19) 
<9td^,Af ^ dtdu weakly in L^/=^(0, T; L^''\Vt- M^)) . (5.20) 

for any z/ G (0, 1/2) , s > 3, where fl535D follows from fICTID . Note that the M- 
projection is kept in the convective term in the limit N ^ oo. 

Applying the Aubin-Lions compactness lemma (cf. [22]), we deduce that, 

0N,Ai Om Strongly in LP((0,T) x fi) (5.21) 

for any p G [1,5/3). Moreover, using fl5.19p . f l5.2Up . a simple interpolation argument, 
and the Aubin-Lions lemma, we obtain that 

V^ti;v,A/ ^ ^xdM strongly in L^((0, T) x fi; M^^^) for G [1, 10/3). (5.22) 

Next, using fl5.12p . f l5.13p . standard interpolation and embedding properties of Sobolev 
spaces, and the Aubin-Lions lemma, we arrive at 

Un,m um strongly in L'((0, T) x Q; R^) , (5.23) 

for some s > 5. Combining this with f l5.22p . we finally obtain 

Un,m ■ VxdN,M Um ■ VAi strongly in L^((0, T) x fi; R^) (5.24) 

for some q > 2. Moreover, from ( I5.13P and ( I5.19p . we have 

dN,M ■ V^.UAr,A,/ -> d^j " V^u^/ Weakly in LP{{0, T) x Q; R^) 
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for some p > 2 , whence 

dtdN,M dtdM weakly in L^{Q, T; L^{Vl- W^)) . (5.25) 
Finally, we have that 

|V,.u^,Mr-2v,.u,v,M ^ \V.umY-^V.um weakly in r'/'-^((0, T) x fi; M^xS)^ 
We conclude that the limit quantities um , , Om, and pm solve the problem 

UM{t, ■) ■ V^<f = for a.a. t e (0, T) (5.26) 



/ 

Jn 



'n 

for any test function if G C°°(f2): 



/ / (wm [itMlM^iiA/ : Va^v?) +PMdiv<^ = / / (§A/ + cr]j/) : V^v^ (5.27) 
Jo Jq ^ ' Jo Jn 

- I uq- ip{0, ■) + T7 / |Va;nM|'^"^V^.WAf : Vr^ip, 
Jn ^^J- Jn 

for any if G Co°°([0, T) xTt-M?), if- n\9^ = 0, where 

(^M = - {^xdM VAi) - {AdM - f{dM)) ® dM ; (5.28) 

and 

§Af = K^Ai) ( ^ — 1 • (5.29) 

Letting -> oo in the equation for d^^M we get 

dtdu + um ■ "^xdM - du ■ V^^um = A^m - /(^m), a.e. in (0, T) x Vt, (5.30) 
supplemented with 

^x{dM)i ■ = 0, z = 1, 2, 3, (5.31) 

dM(0, ■) = do,M. (5.32) 

The passage to the limit in (15. 5 p is more delicate. Actually, the weak lower 
semi- continuity of convex functionals on the right-hand side gives rise to 

dtOM + div {OmUm) + divg^l,^ > ^\V,umY + Sa^ : V^Um + IA^a/ - /(^a/)!' (5.33) 
satisfied in the sense of distributions, with 

qii ■ n\an = 0, (5.34) 

^m(0,-) = ^o,A^, (5.35) 

where 

Qm = -HOm)'^xOm - h{dM)dM{dM ■ ^xOm)- 
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Next, we claim that the total energy is conserved, namely 

dt J^{^\um\^ + ^^l^^M^ + F{dM)) =0. (5.36) 
Indeed, combining (15. ip with v = uj^^m and (l5.5H5.6p . we obtain 



n V2 



J (iyxdN,M ^xdN,M) ■ VxUAr,A/ + {{AdN,M " /(fiiV,A/)) ® C^JV,A/) : V^Uat^m 



+ \AdN,M - f{dN,M) 

whence, by virtue of (15.21) and after a straightforward manipulation, we get 

dt f (^-{un^mI"^ + On,m + -| Vx^at^mI^ + F{dN,M)^ = 0, 

yielding, by passing to the limit as N ^ oo, the desired conclusion (I5.36p . 

Now, we want to show that (I5.33P is actually an equality. Taking v = U7v,a/ in 
(15. ip we get 



(5.37) 



1 1" 



'n 



obtain 



- I|uo|Il2(q) + 2 / / 
Jo Jn 

Next, thanks to (l5.12f[57T4|) . we can take um as a test function in (15.271) to 

+ o/ f^{OM)\^xUM + ViuM\^ (5.38) 
^ Jo Jq 



2 '■^ 



+ I I |V^ltA/|''"^V^UAf : V^UAf 



Jn 

T 



a A/ • ^xUm ■ 



= l|uo||i.(f,) + 2 /" I 
Jo Jn 

Now, multiplying (15. 2 p by AcZat a/ — fidjy^M), we obtain 

||V,d^,A/(t)||i2(f,)+ / F{d^^M){t) + 2 [ [ \AdN^M- f{dN.M)? (5.39) 
Jn Jo Jn 

= ||Vxdo|lL2(f7)+2 / F{dQ)+2 I {mn^m-^ xdN,M — dN,M-^ x'^N,M^ ^dN,M~ f^dNAl))- 

Jn Jo 
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|2 



(5.40) 



Analogously, multiplying fl5.30p by Ad^ — f{dM) we get 

/ F{dM){t) + 2 / / |AdM-/(dM) 

Jn Jo Jn 

= II Vxdo||i2(f^) + 2 [ F{do) + 2 [ (um ■ VAi - d&i ■ V^.ua/, Ad^ - /(dA/))- 
Jn Jo 

Taking the sum of (15.371) and (I5.39P and 05.381) with (15.401) . and, finally, passing to 
the limit as — oo , we obtain 

/ IV^Utv^a/T -> / / |VxUAf|^-2VxUAf : V:,;Ua^, 

Jn Jo Jn 



/ / \AdN,M - f{dN,M)\^ ^ [ [ \AdM-f{d, 

Jo Jn Jo Jn 

entailing, by means of standard Minty's trick and monotonicity argument, 

V,u;v,M ^ V,ua, strongly in U{{Q,T) x fijM^xS), 

AdN,M Adu strongly in L2((0,T) x VL-M?). 
Consequently, the inequality (I5.33P may be replaced by the equality 

^ 1 9 

dtOM + div {OmUm) + div qij = — | V.^A/r + : ^xUm + \AdM - f{dM)\ ■ (5.41) 

Taking ua/V^, with G ^^((0, T) x fi) , as a test function in (I5.27p . testing (I5.30p 
by ^^^ip, adding both relations to (I5.4ip multiplies by (/J, and using (I2.17p . we get 
an M-analogue of (13. 8p . namely: 

dt Q|ua^|^ + eAf^ + div Q|ua/|^[ua/]a/ + ca/Ua/^ (5.42) 

+div [pm^m + Qa/ - Sa/Ua/ - ctm ua/) = div {v ^dM ■ {AdM - f{dM))^, 
with the internal energy 

e-M = -|Va;dA4"|^ + F{dM) + 

and the fiux 

Qm = —k{6M)yx(^M — h{6M){dM ■ yx0M)dM — AVx^A/ ■ V^Ua/ ■ dM- 
Finally, we can multiply (15. 5 p by H'{6m)^, obtaining 

f [ H{dM)dt^+ f I {H{9M)viM + H'{dM)q'it)-V.^ (5.43) 
Jo Jn Jo Jn 
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+ \AdM - /(dA/)r) + H"{9M)qi ■ V.^M^y^ 

- / H{9omMO,-), 
Jn 

for any if G C^([0,T) x > 0, and any smooth, non- decreasing and concave 

function H. To be precise, we have however to remark that, at this level, we do 
not have sufficient regularity in f l5.5p to use H'{6m)^ directly as a test function. 
Nevertheless, the procedure could be justified by a standard regularization argument 
and then taking the (supremum) limit. This is also the reason why we get the < 
sign, rather than the equality, in (15. 43 p . This concludes the passage to the limit for 
N ^ oo. 



5.2 Passage to the limit as M — oo 

Our final goal is to let M ^ oo in and f05]) . We notice that the 

limits in f l5.12p . f l5.16f[5?22|) still hold when letting M — > oo. On the other hand, we 
now have 

dtVLM dtVL weakly in L^(0,T; M^)) , (5.44) 

Pm P weakly in L^((0, T) x 1]) , (5.45) 
dtdu^dtd weakly in L2(0,T;L3/2(f^;M3)), (5.46) 



and, obviously. 



M-i/("-^V^ua/ ^ strongly in ^"-^((0, T) x fi). (5.47) 



The above relations are sufficient to pass to the limit M — )■ oo in (I5.26f[5732|) 
to recover (I3.5H3.7I) . In addition, by (15.111) and the previous estimates, we get 



-+Pm] Um \ bounded in L'((0, T) x fi; R^) for some i > 1 , 
/ J A/>0 

{^MUA/}Af>o bounded in L\{Q,T) x Vt^W?)) for any q G [1, 10/9) , 
{cj]jfuM}j^^g bounded in V{{Q,T) x fi; M^) for some i > 1 , 
qA/ bounded in bounded in -^^''((0, T) x VL; M^) for some l> 1 . 

Notice that we used here in an essential way the fact that rjir— 1) > 10/7. 

As a consequence, we can pass to the limit in (I5.42p and to the lim sup in (I5.43P 
(thanks also to the positivity and convexity of the terms on the last line of fl5.43p ) to 
deduce the desired conclusions (13. 8 p and (13. 9p . This completes the proof of Theorem 

To conclude, we remark that the above estimates are not sufficient for passing 
to the limit in (I5.4ip with respect to M — )■ oo , due to the lack of strong convergences 
of the terms appearing on the right hand side. 
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